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ABSTRACT 


A  class  scheduling  problem  is  formulated  as  a  graph 
coloring  problem.   A  computer  program  based  on  an 
algorithm  recently  developed  by  Welsh  and  Powell,  Computer 
Journal,  Vol.  10,  May  1967,  pp.  85-86,  is  used  to  obtain  a 
solution  to  the  coloring  problem.   While  the  program  fails 
to  provide  an  acceptable  schedule  in  this  application,  the 
results  indicate  that  improvements  in  the  coloring  algorithm 
may  yield  acceptable  schedules. 
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1.   Introduction. 

The  procedure  of  scheduling  classes  varies  with  each 
academic  institution.   In  most  instances  schools  announce 
beforehand  a  list  of  courses  and  times  they  are  to  be 
offered.   In  this  case  students  schedule  themselves  to 
courses  which  do  not  conflict.   A  different  scheduling 
problem  occurs  in  the  case  where  students  are  required  to 
take  .predetermined  courses.   Here  the  problem  is  to 
schedule  the  courses  so  that  students  have  no  conflicts. 
The  problem  of  scheduling  classes  at  the  Naval  Postgraduate 
School,  Monterey,  California,  falls  in  this  category.   It 
is  this  problem  which  is  investigated  here. 

The  digital  computer  has  been  used  as  an  aid  in 
scheduling  before.   Appleby,  Blake,  and  Newman  ["lj  attempt 
to  produce  school  time-tables  using  a  general-purpose 
digital  computer.   Their  suggested  methods  of  approach 
include 

(1)  a  random  trial  solution  updated  to  eliminate 
conflicts, 

(2)  a  trial  of  all  combinations  possible  until  one 
with  no  conflicts  is  found, 

(3)  a  random  buildup  of  entries  into  a  schedule, 
adding  entries  only  when  conflicts  do  not  occur,  modifying 
previous  entries  to  make  new  entries  possible,  and 

(4)  a  heuristic  approach. 

The  importance  of  knowing  whether  or  not  a  solution  does 
exist  is  stressed.   The  last  of  three  successively  more 
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complex  programs  is  estimated  to  speed  up  scheduling  by  a 
factor  of  200  to  400  over  the  hand  operation. 

Csima  and  Gottlieb  £6J   describe  a  time-table  construc- 
tion technique  which  uses  matrices  of  zeroes  and  ones  to 
produce  results  in  the  limited  situations  attempted.   In 
real  situations  where  preassignments  or  special  constraints 
appear,  further  progress  must  be  made.   It  is  also  hoped 
that  with  refinements  in  the  program  the  computer  execution 
time  can  be  reduced. 

A  method  of  determining  examination  time-tables  has 
been  described  by  Cole  £5j.   Citing  Sherman's  Q9]  pessimism 
in  finding  a  schedule  for  a  large  school,  he  emphasizes 
producing  an  acceptable  solution  for  as  many  courses  as 
possible  with  reasonable  computer  time.   Dependent  on  the 
imposed  constraints,  the  execution  time  on  the  small-store 
computer  was  estimated  to  be  several  hours  if  all  of  340 
courses  were  considered. 

Another  method  using  GASP,  Generalized  Academic 
Simulation  Program,  a  program  written  at  M.I.T.  ^1~], 
approaches  the  problem  by  assigning  students  to  a  previously 
determined  course  schedule. 

This  paper  will  treat  the  scheduling  problem  by 
assigning  days  and  hours  to  courses,  so  that  the 
constraints  imposed  by  the  students  as  well  as  instructors 
are  satisfied. 
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At  the  Naval  Postgraduate  School  there  are  aspects  to 
the  problem  which  make  it  different  from  most  institutions. 
These  aspects  are  treated  as  restrictions  to  the  basic 
scheduling  problem  and  will  have  an  effect  on  its  solution. 
For  instance,  the  government  allots  a  given  length  of  time 
to  students  to  follow  a  prescribed  curriculum  of  study. 
Since  certain  courses  are  taught  in  sequence,  students 
are  required  to  take  certain  courses  at  certain  times. 
The  students'  schedule  thus  loses  some  flexibility. 
Another  aspect  to  be  considered  is  the  fact  that  all  naval 
aviators  must  be  scheduled  either  a  complete  morning  or 
afternoon  for  weekly  aircraft  flights.   Other  restrictions 
include  days  off  for  instructor  research  and  allowances 
made  for  faculty  and  staff  education  at  the  school  itself. 

A  graph  theoretic  approach,  which  makes  use  of  the 
computer  for  the  necessary  calculations,  is  used  to  attempt 
the  solution  to  the  scheduling  problem.   A  graph  theoretic 
model  of  the  scheduling  problem  is  constructed.   The 
problem  of  scheduling  courses  without  conflicts  is  shown  to 
be  equivalent  to  a  graph  coloring  problem.   A  digital 
computer  program  is  produced  to  solve  the  coloring  problem. 

Section  2  of  this  paper  defines  the  problem  investi- 
gated.  Preliminary  definitions  and  terminology  are  first 
presented  in  graph  theoretic  terms.   The  graph  theoretic 
model  is  then  constructed  with  these  tools.   In  section 
3  an  algorithm  for  graph  coloring  is  analyzed.   This 
algorithm  then  is  used  to  construct  a  coloring  program. 


The  results  of  this  program  are  displayed  in  Table  1. 
Observations  and  conclusions  based  on  these  results  are 
found  in  Section  4.   The  Appendix  contains  both  general 
and  detailed  flowcharts,,  which  explain  the  coloring 
program  used. 
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2.   Definition  of  the  Problem. 

A.   Preliminary  Definitions. 

We  wish  to  consider  the  basic  problem  of  scheduling 
from  a  graph  theoretic  standpoint.   But  first  we  need  the 
following  basic  definitions  taken  from  Berge  £73j^  Chapters 
1  and  2. 

Definition  1.   Given  the  sets  X  and  Y.   A  function  f  , 
mapping  X  into  Y,  is  multi -valued,  if  to  each  xeX  there 
is  a  well  defined  subset  y  x^Y.   Note:  ]/   x  may  be  the 
null  set. 

Definition  2_.  yx  is  said  to  be  the  set  of  followers  of 
x. 

Definition  3.   A  pair  denoted  by  G=(X,r  )  is  said  to  be  a 
graph  whenever 

(l)   X  is  a  set, 
and   (2)  ^  is  a  multi-valued  function  mapping  X  into  X. 
It  is  customary  to  represent  the  elements  of  the  set  X  as 
points  in  a  plane. 

Definition  4.   If  G=(X,r  )  is  a  graph  and  x6  X,  then  x  is 
said  to  be  a  point ,  vertex,  or  node  of  the  graph  G. 
Definition  5_.   If  x€  X  and  y€  X  such  that  y£  f   x,  the 
ordered  pair  (x,y)  is  said  to  be  an  arc  of  the  graph  G. 
Arcs  of  a  graph  are  usually  represented  by  an  arrow  from 
the  first  to  the  second  element  of  the  ordered  pair. 

In  case  the  order  in  which  the  vertices  appear  in  the 

ordered  pair  (x,y)  has  no  significance,  we  denote  the  pair 

by  Cx>yJ  =  [j  >^Q  •   Such  an  unordered  pair  is  called  an  edge 

of  the  graph  G. 
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An  edge  of  the  graph  is  usually  represented  by  a 
curve  joining  two  vertices. 

Definition  6.   If  G=(X,f)    is  a  graph  and  U=/[x,y]  :y6  iTx 
for  all  x  j   the  set  of  edges  of  G,  then  an  alternate 
notation  for  the  graph  is  the  pair  (X,U).   Since  this 
graph  contains  only  undirected  edges,  it  is  sometimes 
referred  to  as  an  undirected  graph. 

Associated  with  every  graph  is  a  matrix.   We  shall 
define  it  as  follows: 

Definition  7.   If  G=(X,U)  is  a  graph,  then  A=(a,.)  is 
called  the  matrix  associated  with  G  if  a.  .  equals  the 
number  of  edges  from  node  x.  to  node  x.. 
Note:   A  is  a  square  symmetric  matrix. 
Definition  8.   Given  the  graph  (X,U),  where  x€  X  and 
y6X,  x^y,  x  and  y  are  said  to  be  adjacent  if  fx,y]  = 

Definition  9.  If  G=(X,U)  is  a  graph  and  if  x6X,  the 
number  of  edges  having  an  endpoint  at  x  is  said  to  be 
the  degree  of  x. 
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B.   Definition  of  the  Scheduling  Problem. 
Definition  10.   A  section  is  a  group  of  students  all  of 
whom  take  the  same  courses. 

Sections  are  used  to  simplify  the  scheduling  problem. 
Determining  a  single  section  schedule  is  easier  than 
finding  a  separate  schedule  for  each  individual. 

Denote  sections  as  s. ,    i=l,..., n, ,  courses  as  c,, 

J-  J.  J 

j=l,...,np*  and  professors  as  p,  ,  k=l,  .  .  .  ,n~..      We  shall 
consider  the  set  X=Js.  i  {j\c  •  \  U  "Ip-  (  as  a  set  of  nodes  for 
a  graph.   We  shall  also  consider  the  sets  of  all  pairs 
U,  =  <fs.  ,c  ."] :  section  s.  must  take  course  c  .  (  and  Uo^-ZTc., 
p,  ^J :  instructor  p,  is  assigned  to  teach  course  c.f. 
Define  U=U, U  Up .   We  see  that  U  with  the  set  X  defines 
the  graph  (X,U). 

Definition  11.   The  grapn  G=(,X,u)  described  above  is 
called  the  data  graph. 

We  shall  consider  an  example  with  the  following 
conditions : 

(1)  Section  s,  is  required  to  take  courses  e,,Cp, 
and  c». 

(2)  Section  s?  is  required  to  take  courses  c.-,  and 

(3)  Instructor  p,  is  assigned  to  teach  courses 
c,  and  Co. 

(4)  Instructor  p?  is  assigned  to  teach  courses 
Cp  and  c^,. 
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Then  as  defined  above,  X=js1, s2 ,c1 ,Cg ,c«, c^,p1,p2Vand 
u"/CBi'°l3'  &1,c2],[s1,c3],fs2,cJ,  Cs2'c3]'Ul'cJ' 

Epl'c3]'[p2'c2]'fp2'c4j/-   The  data  graph  G=(X'U)  ls 
illustrated  as  follows  in  Figure  1. 


p* 


Fig.  1.  Example  Data  Graph 
Definition  12 .   Given  the  data  graph  G=(X,U)  described 
above,  a  conflict  is  said  to  occur  between  classes  c. 
and  Qy    i^j, 

if  (l)  [s,  ,c,"]  €  I J  and  fs  ,  cJ^U  for  some  k, 
or  if  (2)  [c^p-jGtJ  and  [c.jpJ^U  for  some  1. 
This  is  the  same  as  saying  these  classes  c.  and  c, 
cannot  be  scheduled  simultaneously  because: 

(1)  the  same  section  s.  must  take  both  classes  and 
could  not  be  present  at  two  places  simultaneously,  or 

(2)  the  same  instructor  p,  must  teach  both  classes 
and  could  not  do  so  simultaneously. 
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We  consider  now  the  concept  of  a  conflict  graph 
deduced  from  the  data  graph. 

Definition  13.  Given  the  data  graph  (X,U),  a  new  graph 
(C,U!)  is  called  a  conflict  graph  if  C=4c.  :c.  is  a  node 
corresponding  to  class  c.  for  all  i  V  and  U'=?<  fc.  ,c  H : 


classes  c.  and  c.  have  a  conflict 


)• 


Note:  c.  cannot  have 


a 


conflict  with  itself;  therefore,  for  all  I,  [c.  ,  c-"]^  U  ' 

The  associated  matrix  with  the  conflict  graph  we 

shall  call  the  conflict  matrix.   The  conflict  matrix 

M=(m. . ),  i, j-1, „ . .,npj  is  such  that 

m.  .=  Co,    if  c.  and  c.  have  no  conflict. 
J  j  -^      J 

/  1.  if  c.  and  c.  have  a  conflict. 
If  we  consider  the  same  example  as  before,  we  can 
display  the  conflict  graph  in  Figure  2  below. 


Fig.  2,  Example  Conflict  Graph 
The  conflict  matrix  in  this  case  is  M=/£)/  /  Q 

/0/  I 
/  /  01 
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Coloring  a  graph  is  a  process  of  assigning  colors  to 
each  node. 

Definition  14.   Given  an  integer  p,  a  graph  G  is  said  to 
he  p -chromatic ,  if  the  nodes  can  be  colored  with  p 
distinct  colors  so  that  no  two  adjacent  nodes  are  assigned 
the  same  color. 

Definition  15.   The  smallest  number  p  for  which  a  graph 
is  p-chromatic  is  called  the  chromatic  number  of  the 
graph  G  and  is  written  k(G). 

Given  a  group  of  classes  with  arbitrary  conflicts 
among  them,  we  solve  the  problem  of  scheduling  by 
assigning  time  periods  to  these  classes.  Where  a 
conflict  arises  between  two  classes  we  require  that  they 
be  assigned  different  time  periods.   On  the  other  hand, 
if  we  consider  these  classes  as  nodes  of  a  graph  and 
their  conflicts  as  edges,  the  problem  of  coloring  the 
nodes  of  this  graph  is  seen  to  be  equivalent  to  the 
above  assignment  of  time  periods. 

The  problem  posed  here  is  not  necessarily  one  of 
determining  the  chromatic  number  of  a  conflict  graph. 
Our  problem  is  to  find  a  suitable  coloring  of  the  conflict 
graph,  which  is  consistent  with  imposed  restrictions 
which  are  now  discussed. 

The  Naval  Postgraduate  School  uses  the  quarter  system 
for  teaching  courses.   Each  course  lasts  twelve  weeks. 
During  the  twelfth  week  of  each  quarter  final  exams  are 
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held.   The  first  eleven  weeks  have  the  same  schedule. 
Normally  each  course  will  meet  several  hours  each  week 
depending  upon  the  nature  of  the  course. 

Earlier  we  assigned  a  node  c.  to  each  course.   We 
shall  now  introduce  nodes  with  two  subscripts  to 
differentiate  between  the  hours  of  a  given  course. 
Consider  a  three  hour  course  c.„   Then  the  hours  or 
nodes  associated  with  these  hours  are  denoted  c.  ,, 
c.  0,    c.  oj  or  simply  c.  .,  j=I,2,3»   If  a  course  is 
a  one  hour  course,  we  shall  let  j=0. 

In  some  cases,  if  there  is  a  large  number  of 
students  taking  the  same  course  in  a  given  quarter, 
it  may  be  necessary  to  divide  this  group  into  two  or 
more  subgroups  for  teaching  purposes.   These  subgroups 
are  assigned  to  take  different  sessions  of  the  same 
course.   These  different  sessions  we  shall  call  segments 
of  that  course.   We  shall  introduce  nodes  with  three 
subscripts  to  differentiate  between  segments  of  a  course. 
In  the  example  above,  if  course  c.  has  two  segments,  then 
the  hourly  nodes  associated  with  the  course  are  c.  .  ,  , 
i  denoting  the  course,  j  being  the  hour,  and  k  identifying 
the  segment.   In  this  case  they  are  c.  .  .  ,  j=l,2,3^ 
k=l,2.   If  a  course  only  has  one  segment  we  shall  let 
k=0. 

Now  we  let  X=js.j(j(c  .  ^  ^  1  (J  fpl    for  all  i,  j,  k, 

1,  and  m.   We  also  consider  U,^/rs.,c.  ,  /7  :  section  s. 

.1  |L  i   ,]  ,k,  U  i 

must  attend  hour  c.  ,  .,  f  and  U^-lYc  .  .  ,,p  1:  instructor 

j,k,lj      2  j1-  j,k,l'tmJ 
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•p   is  assigned  to  teach  hour  c.  ,  ,  f.   If  U=U,L/lU,  then 
we  have  G-=(X,U)  as  our  data  graph.   We  assume  that  as  the 
need  for  segments  arises,  the  appropriate  section  to 
segment  and  segment  to  professor  assignments  are  made 
by  the  academic  departments.   We  may  thereafter  regard 
these  segments  as  different  courses. 

Given  our  new  data  graph,  we  define  conflicts  in 
exactly  the  same  manner  as  In  Definition  12,  except 
that  we  now  allow  course  nodes  to  have  three  subscripts. 

We  now  consider  the  previous  example  with  the 
modifications  below. 

(1)  Let  s,  take  courses  c,,  c?,  and  the  first 
segment  of  course  o_. 

(2)  Let  Sp  take  course  c^  and  the  second  segment 
of  course  c^. 

(3)  Courses  c,j  Cg,  and  c~  will  all  be  three  hour 
recitation  courses. 

(4)  Course  c^  will  consist  of  three  recitation 
hours  and  two  laboratory  hours  to  be  scheduled  together. 

(■■))   Instructor  p..  will  teach  both  segments  of  z    . 
In  this  case,  5C--fc1)1)0,c2jJj0,c,iJ;J,c4ik;0-  1=1,2,3, 
j=l,2 ,k.--=l, . . .  ,5  l.      Figure  3  illustrates  G,  the  data 
graph. 

Department  meetings,  instructor  research  time, 
aviator  flight,  time,  and  other  sued  periods  imposing 
constraints  are  dii     ed  on  the  conflict  graphs,  in 
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Pig.   3.  Data  Graoh 
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the  case  of  a  one  hour  lecture  for  all  students,  we  add 
to  X  a  lecture  node  c,  Q  Q  for  some  1,  and  to  U  we  add 
edges  £c,  n  r)>^~]   for  aH  *•   In  tne  case  of  a  one  hour 
meeting  for  all  professors  in  a  given  department,  we 
add  to  X  a  node  c,  n  n  for  some  k,  and  to  U  we  add  the 
set  of  edges  I  £c,  Q  Q,p  J  :  instructor  p.  belongs  to 
that  departments   Similarly  other  constraints  are 
represented. 

Our  data  graph  is  G=(X,U).   X  is  the  set  of  all 
nodes  denoting  sections,  professors,  courses  and  their 
segments,  flights,  instructor  research,  and  special  hours 
such  as  department  meetings  or  special  lectures.   The 
set  U  consists  of  all  applicable  edges. 

Using  the  above  data  we  must  introduce  a  slightly- 
different  concept. 

Definition  16.   Given  the  data  graph  G=(X,U),  a  new  graph 
G'=(C,U')  is  called  an  hourly  conflict  graph  for  G,  if 

C  =  /c.  .  .  :  c.  .  ,  is  a  node  corresponding  to  course  c, 
(_J-jJ>K--LjJjK-  1 

hour  j,  and  segment  k  (  and  U'  =  <[c.  ,  ,  ,c,    "):  hours 

J        J_uiJj,k   l,m,nJ 

c.  .  .  and  c,     have  a  conflict] . 
i,j,k      l,m,n  J 

Based  on  our  previous  example  the  hourly  conflict 
graph  is  illustrated  in  Figure  4. 
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Fig.  4.  Hourly  Conflict  Graph 

Given  a  coarse  that  meets  several  times  a  week,  we 

consider  now  conflicts  between  the  hours  of  the  course 

described  above.   These  hours  may  be  divided  Into 

recitations  or  laboratory  work.   For  a  given  course, 

using  the  second  subscript,  we  shall  number  the  recitation 

hours  first.   For  instance,  a  non-segmented  course  c, 

a  five  hour  course  consisting  of  three  recitations,  will 

have  its  hours  denoted  c.  ,  ^,  c.  „  .,,  c.  n  A.  c.  ,,  ~, 

1,1,0'      1,2,0'      1,3,0'       i,4,0' 

and   c.    ,_    n,    the   first    three   of  which   correspond   to 
recitations.      Note:      Fc.     .   ^,c.    ,     rTl€  IP    for  all   j^k. 
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We  shall  restrict  our  problem  somewhat  by  allowing 

only  one  recitation  per  course  to  be  held  daily.   This 

brings  about  the  following  definition: 

Definition  17.   If  C  is  the  set  of  all  course  hour  nodes, 

then  there  exists  a  daily  conflict  between  c^    .  k  and 

c,     provided: 
l,m,n  ^ 

(1)  1-1, 

(2)  k-n, 

(3)  #m, 

and   (4)   j  and  m  both  correspond  to  recitation  hours. 

We  now  augment  our  information  regarding  conflicts 

by  the  following: 

Definition  18.   If  G=(X,U)  is  the  data  graph,  then 

G' '=(C,U' ' )  is  called  the  daily  conflict  graph  for  G,  if 

C  is  the  set  of  all  course  hour  nodes  and  U'^/fc.  .  .  , 

(ui,  j,kJ 

c.    ,  1:  e,  ,  ,  and  c.    .  have  a  daily  conflict^. 
i,m,kJ   i,j,k      i,m,k  J  J 

If  we  consider  our  previous  example  again,  we  see 
that  Figure  5  below  illustrates  G'',  the  daily  conflict 
graph. 


%o 


cv,o   G/V>   Q)J/0  Q/V>  %  %\    ty,i  G,3,i    %o  %o     Q^ 

Fig.    5.    Daily  Conflict    Graph 
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In.  the  case  of  laboratory  hours  of  a  course  or  flight 

hours  for  aviators,  these  hours  must  be  scheduled  together 

during  a  given  day.   To  consider  this  we  introduce  the 

following: 

Definition  19.   If  G=(X,U)  is  the  data  graph,  then  S= 

(C,tP ' ' )  is  called  the  sequence  graph  for  G,  if  C  is  the 

set  of  all  course  hour  nodes  and  U',,  =  '{rc.  .  ,  ,c,    "7: 

i*-i,J,k*    l,m,nJ 

c.     .   ,     must  be   scheduled   together  with  c,  v. 

The   sequence   graph   S  for   our  previous    example   is 
shown   in  Figure  6  below. 

Gfto  CVo  Cw  Cm  Gf.fii0  CfcC 

»  •  •  #  0 


C'/',*  %*   <V,»   ^s,o  %i   fy  CV)l  QM   fy/0  ^     ^ 

Pig.  6.  Sequence  Graph 
We  have  a  data  graph  which  summarizes  our  input 
information  containing  all  section,  course,  and  instructor 
assignments.   The  hourly  conflict  graph  G'  displays  course 
hours  which  must  be  scheduled  on  different  hours.   The 
daily  conflict  graph  G''  displays  course  hours  which  must 
be  scheduled  on  different  days*   Our  sequence  graph  S 
indicates  which  course  hours  must  be  scheduled  together 
in  block  form.   The  problem  is  to  schedule  these  course 
hours  in  a  five  day  week.   The  school  allows  nine  class 
hours  per  day. 
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We  will  now  consider  the  concept  of  a  color  vector 
with  two  components,  the  first  of  which  will  be  a  color 
corresponding  to  the  day.   The  second  will  correspond  to 
the  hour  within  the  day.   Our  problem  is  to  color  the 
course  nodes  of  the  hourly  conflict  graph  G'  so  that  no 
two  adjacent  course  nodes  have  the  same  color  vector. 
The  coloring  of  the  nodes  of  G'  may  be  thought  of  as  a 
simultaneous  coloring  of  all  three  graphs  G',  G1',  and 
S  subject  to  the  following  restrictions: 

(l)   no  two  adjacent  nodes  of  G'  are  assigned  the 
same  color  vector, 

and   (2)   no  two  adjacent  nodes  of  G"  are  assigned  the 
same  day  color, 

and   (3)   adjacent  nodes  of  S  are  assigned  the  same  day 
colors  and  sequential  hour  colors. 

Note:   It  is  sufficient  for  two  color  vectors  to  be 
different  when  either  or  both  of  the  corresponding 
components  are  different.   In  view  of  this  it  is  interesting 
to  note  that  satisfying  (2)  above  automatically  satisfies 
(l)  for  the  same  two  nodes. 

In  this  coloring  process  we  are  restricted  to  five 
day  colors  since  the  schedule  desired  is  a  weekly 
schedule.   In  addition,  the  number  of  hour  colors  within 
a  day  must  be  restricted  to  nine.   Given  a  satisfactory 
coloring  of  the  conflict  graph,  we  assign  days  to  the 
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day  colors  and  hours  to  the  hour  colors.   The  resulting 
schedule  meets  all  the  constraints  Imposed  provided 
we  have  not  exceeded  the  maximum  number  of  colors 
allowable  in  each  case. 
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3,   The  Coloring  Process. 

A.  The  Algorithm, 

Although  we  cannot  be  certain  of  obtaining  the 
chromatic  number  for  a  given  arbitrary  graph,  coloration 
solutions  may  be  attained  by  trial  and  error  methods  or  by 
linear  programming  techniques  £ 3 J.   An  algorithm  for 
coloring  a  graph,  introduced  by  Welsh  and  Powell  Cioj,  has 
been  used  to  attempt  to  solve  the  vector  color  problem  of 
the  Naval  Postgraduate  School  course  graph.   This 
algorithm  is  described  as  follows: 

Let  G=(X,U)  be  an  arbitrary  graph  where  X  is  the  set 
of  nodes  and  U  the  set  of  edges.   Denote  the  vertices  X= 
«A^ ,    i=l,...,n£.   Let  the  degree  of  the  vertex  A   be 
denoted  by  cL,  where  1=1,.,., n.   Assume  that  the  vertices 
have  been  numbered  so  that  their  degrees  form  a  non- 
increasing  sequence  d^  d^  .  .  .  ^n*      We  now  define  T.CX, 
where  1=1,..., t.   T   is  defined  as  follows: 

(1)   A1€T1, 


and      (2)      if     ^A     }    £    '1!        where   i.-l   and   i.<  i0<  .  .  ,<i 

c       k*-i  L  id 

for    some  m,    then  A.€T.  ,    if   and    only    if, 

(a)      j  >i    , 

and       (b)       [  A    ,A    /]  d    U   for   k-l,...,m. 
j        k 
Note:      T2^0. 

Similarly  we  define  T    ,   where  Tp  C  (X-T    ).      T      is 

constructed  as  follows: 

(l)     If   I   Is   the   least    Integer  for  which  A.^T,, 

then  A.€  T„« 
l        2 


m 
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(2)      If     (a      J6T,   where   J   <  j    <  %  .  ,<i    ,    then 
A   €  T2,    if   and   only   if, 

(a)        />Jp> 
and      (b  )      pA|    ,A.     ]*   U  where  k=l,.,.,p, 

If  Tp=0,    we  are   finished.      If   not,    we   construct   T.C(X- 

Z 

{JT .  )  in  the  same  way  as  we  constructed  Tp . 

Since  X  is  finite,  there  will  exist  a  least  integer 

r-/ 
r  such  that  T=0,      We  note  also  that  (JT.=X  and  that  if 

i^j,  then  5?.  f)T  .=$*      We  can  now  assign  a  different  color 

to  each  T.  .   If  A.^A.^T,,  i^j,  for  some  k,  then  Qa..,A."] 

«U.   Thus  no  adjacent  nodes  will  have  the  same  color. 

The  algorithm  furnishes  an  upper  bound  for  the  chromatic 

number  k(G).   Thus  k(G)^r-l. 

Welsh  and  Powell  further  state  that 

k(G)^  r-l^max  min(i,d.+l) 
Proof  of  this  upper  bound  comes  from  analysis  of  the 
algorithm.   By  construction,  each  node  belonging  to 
(X-T1 )  is  adjacent  to  a  node  of  T  m     By  construction 

Ml 

we  see  that  each  node  of  (X-UT  •)  is  adjacent  to  at  least 

one  node  of  each  of  j  T.(   .   Hence,  A,  6(X-L/T.  )  implies 

L      JJ.  K       j   1 

d,  ^m.  , 

Consider  (X-^/T.).   A.  cannot  belong  to  this  set 
because  d.>d.+l.   Hence, 


Ad.+l, 

(i)  A.ei/Tj. 


Since   each  T.    contains   at    least    one    element, 

J      » 

(2)      A6L/T, 
/       J 
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Together  (l)  and  (2)  above  imply  that 

w) 

(3)   A.6L/T.,  where  0?(i  )=min(i,d  +1 ) . 

1   »  J  <#•/   #*; 

Let  d  (G)=niax  (p(i).   Then  A,  €  L/T  C  L^T  for  all  i. 
I      J  i   j   J   \      J 

Setting  q=  J*(g)+1  we  see  that  T=0.   Thus  for  the  least 
integer  r  where  T=0,  this  algorithm  ensures  r^q.   Hence, 

k(G)^r-l^q-l=  cf (G)=m^x  min(i,di+l).   The  desired 
inequality  is  proved. 

The  above  O  (g)  provides  us  with  an  upper  bound  for 
the  chromatic  number  of  a  graph.  Berge's  ["2J  statement  of 
Brooks'  \Ja\   theorem  on  coloring  the  nodes  of  a  network 
gives  us  insight  about  a  lower  bound.   This  statement, 
simplified  to  apply  to  our  problem,  is  as  follows: 
Theorem  (Brooks ) .   A  graph  G  of  maximum  degree  q  where 
q>2  can  be  colored  with  q  colors  except  in  the  case  where 
a  component  of  G  consists  of  q-f-1  vertices,  each  of  which 
is  joined  to  the  other  q  vertices. 

The  above  mentioned  component  is  referred  to  as  a 
complete  subgraph  of  Q.   In  other  words,  if  G  contains  a 
complete  subgraph  on  m   vertices,  then  k(G)^.rn.   Thus  the 
largest  complete  subgraph  of  G  sets  a  lower  bound  on  k(G). 

To  illustrate  the  use  of  this  algorithm  in  coloring, 

consider  the  given  arbitrary  graph  in  Figure  7  below. 

I 


Fig.  7.  An  Arbitrary  Graph 
28 


Numbering  of  the  nodes  must  be  done  so  that  their 
degrees  form  a  nonincreasing  sequence.   One  such 
renumbering  is  listed  below. 


Node 

Degree 

R< 

^numbered 

1 

3 

A3 

2 

4 

Al 

3 

3 

A4 

4 

4 

A2 

5 

3 

S 

6 

3 

A* 

Figure  8  below  illustrates  the  arbitrary  graph  with 
nodes  relabeled.  hi 


Fig.  8.  Relabeled  Arbitrary  Graph 
Using  the  algorithm  to  form  the  disjoint  sets  T., 
i=l„j4jWe  see  the  following: 

T1={A1'A2^  T2  =  fA3'A5}>  T3^4^A6^  and  T4=ef* 

To  obtain  a  satisfactory  coloring,  we  assign  the  nodes  of 

T,,  Tp,  and  T^  colors  1,   2,   and  3  respectively. 

We  note  from  the  above  that  this  arbitrary  graph  is 

3  colorable  by  the  algorithm.   Of  interest  is  a  comparison 

of  this  number  with  the  lower  and  upper  bounds  of  the 

chromatic  number  for  this  graph.   By  observation,  the 
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largest  complete  subgraph  has  3  nodes;  therefore,  the  lower 
bound  is  3.   Since  the  number  of  colors  in  our  solution 
equals  the  lower  bound,  we  deduce  that  the  chromatic 
number  is  in  fact  3.   The  upper  bound  of  k(G),  cT(G)=max 
min(i,di+l)=maxjl,2, 3,4,4,4 f is  4  for  this  example. 

While  in  general  we  will  not  know  what  k(G)  is  for  an 
arbitrary  graph,  this  algorithm  does  ensure  a  p-coloring 
where  p^  J"(G)=max  min  (i,d.+l). 
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B.   Coloring  the  Graph  by  Color  Vectors. 

In  coloring  our  conflict  and  sequence  graphs  we  are 
restricted  to  five  day  colors  and  nine  hour  colors  per 
day.   Flowchart  A  in  the  Appendix  illustrates  the  coloring 
procedure  and  application  of  the  Welsh  and  Powell  algorithm 
to  this  problem  by  means  of  a  digital  computer  program 
using  Fortran  IV  language.   This  procedure  is  described 
briefly  in  chronological  order  as  follows: 

(1)  Section  and  instructor  course  assignments  are 
read  in, 

(2)  A  course  directory  is  established. 

(3)  An  indexed  list  of  conflicts  is  determined  based 
on  section  and  instructor  conflicts. 

(4)  Course  hour  data  is  read  in, 

(5)  Course  hour  nodes  are  given  day  colors  using  an 
arbitrary  selection  process. 

(6)  Course  hour  nodes  for  each  day  are  given  hour 
colors  using  the  algorithm. 

Data  for  the  program  has  three  forms.   First  are  the 
sections  cards  containing  the  section,  the  number  of 
students,  the  number  of  aviators  for  that  section,  and 
each  course  taken  by  that  section.   To  differentiate 
between  courses,  segmentation  and  lab  designations  were 
used.   To  indicate  instructor  assignments,  data  cards  were 
used  with  the  same  format  as  on  section  course  assignment 
cards.   This  procedure  simplified  the  computation  of 
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conflicts  between  course  hour  nodes.   Flowchart  B  in  the 
Appendix  illustrates  the  read-in  and  print-out  of  conflict 
data. 

Flowchart  C  in  the  Appendix  illustrates  the 
computation  of  a  course  directory.   Each  course  is 
associated  with  an  integer,  and  the  integer  listing  of 
these  courses  is  used  for  all  further  computations.   This 
course  listing  is  printed  out  for  reference.   Included  in 
this  computation  is  a  random  number  generator  used  for 
random  reordering  of  this  course  directory.   Flowchart 
I  in  the  Appendix  illustrates  RAN  the  generator  function 
used. 

The  list  of  conflicts  is  computed  next.   This  is 
shown  in  Flowchart  D  in  the  Appendix.   In  the  example 
described  in  Section  2,    the  directory  assumes  the 
following  form. 

cc(i)      I 

C0URS1         1 

C0URS2         2 

C0URS31       3 

C0URS32       4 

C0URS4        5 

C0URS4L       6 
C0URS31  and  COURS32  are  segments  1  and  2  respectively  of 
C0URS3  and  C0URS4L  is  the  lab  portion  of  C0URS4.   The 
conflict  list  for  the  above  example  is  shown  below. 
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J  1234567 

INDEX(j)         1      4      8   11   15    18  21 

I  1     2     3     4     5     6     7     8     9  10  11  12   13  14  15    ... 

LIST  (I)  234135612413562... 

The  above   Information   is   interpreted   as   follows: 
C0URS1  whose   integer   representation   is   1   has   conflicts 
with  C0URS2,    C0URS31,    and   C0URS32,    whose   integer 
representations  are  2,    3,    and   4   respectively.      The 
conflicts   are  found  by   entering  the   list   at   index(l)   and 
reading  all   entries   up   to   index(2),   where  the   conflicting 
courses    of   C0URS2   begin.      This   notation  has   the  advantage 
of   requiring   less  memory   space  and   less  print -out   space 
than  the   comparable  matrix  representation  discussed 
earlier.      The   nonzero   elements    of   the  I-th  row    or 
column   of   this  matrix  are   listed   starting  at   index(l). 

The   list  which  describes   the   conflicts  between  the 
nodes   is   ordered   so  that   the   nonzero   elements   corresponding 
to   the   I-th  row    or   column   of   the  matrix  form  an   increasing 
sequence.      Flowchart   E   in  the  Appendix  describes   this 
procedure  and   the  print-out    of   its   results.      Flowchart 
J   in  the  Appendix  displays   the   subroutine  RELIST   used   in 
this  procedure. 

Incidental   to  the  main   objective   of   the  program,   an 
upper  bound   for   the   chromatic    number   of   the   course  graph 
is   computed.      During   the  development    of   the  program  it 
served   to  give  an  indication   of   the  degree   of   difficulty 
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to  be  encountered  later  in  the  coloring  phase.   Flowchart 
F  in  the  Appendix  shows  the  computation  of  the  upper  bound 
for  this  graph  using  courses  as  nodes. 

We  next  read  course  hour  information  into  our  program, 
in  order  that  day  colors  may  be  assigned  to  the  course 
hour  nodes.   Flowchart  G  in  the  Appendix  illustrates  this 
input  and  assignment  of  day  colors.   Flowchart  K  describes 
the  subroutine  DAYCOL  called  for  assignment  of  day  colors. 
This  routine  assigns  recitation  hours  of  the  same  course 
different  day  colors.   Each  lab  course  is  recognized  as 
such  in  order  that  its  hours  may  be  scheduled  during  the 
same  day.   A  listing  with  index  is  compiled  as  shown,  to 
remember  these  day  color  assignments.   In  our  example  the 
course  hour  information  is  as  follows: 
I      CC(I)        HOURS 

1  C0URS1     3  recitation 

2  C0URS2     3  recitation 

3  C0URS31    3  recitation 

4  COURS32     3  recitation 

5  C0URS4     3  recitation 

6  C0URS4L    2  lab 

Using  the  above  course  hour  information  DAYCOL 
assigned  the  day  colors  in  our  example  as  follows: 
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Day  Color 

1 

2 

3 

4 

— — T 

5 

Course 

1 

1 

1 

2 

2 

Hour 

2 

3 

3 

3 

4 

Nodes 

4 
6 
6 

4 

5 

5 

5 

To  each  day  we  now  apply  the  coloring  algorithm. 
Flowchart  H  in  the  Appendix  gives  the  general  procedure 
for  the  hour  coloring.  Flowchart  L  illustrates  the 
HRCOLR  subroutine  used  in  that  procedure.  When  labs 
having  more  than  one  period  are  scheduled  on  any  given 
day,  the  integer  representation  for  that  course  will 
be  repeated  for  that  day.   To  enable  these  nodes  to 
maintain  their  separate  identities,  a  renumbering  of 
nodes  is  made.   Relative  to  this  new  numbering,  an  indexed 
list  of  conflicts  is  computed  from  the  master  list  for 
each  particular  day.   Then  the  degree  of  each  node  is 
determined.   Knowing  this,  the  program  assigns  colors  to 
each  node  according  to  the  coloring  procedure  described 
earlier.   The  result  is  the  second  or  hour  component  of 
a  color  vector.   Having  day  colors  already  determined  by 
DAYCOL,  the  process  is  complete. 

Using  the  above  program,  data  from  the  second 
quarter  of  the  I967-I968  academic  year  at  the  Naval  Post- 
graduate School  was  processed.   No  special  periods  such  as 
flight  hours,  instructor  research  time,  or  departmental 
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meetings  were  included  in  this  data.   Included  were  424 
courses,  403  sections,  and  194  instructors.   The  number 
of  conflicts  for  any  course  varied  from  0  to  52;  however, 
over  half  of  the  courses  had  conflicts  ranging  from  5  to 

13. 

The  number  of  hour  colors  required  for  the  five  days 
were  9,  8,  10,  10,  and  10.   To  investigate  the  effect  of 
the  DAYCOL  assignments  on  the  total  number  of  hour  colors, 
a  random  number  generator  was  used  to  reorder  the 
directory  of  courses  prior  to  the  DAYCOL  routine  in  the 
program.   Flowchart  J  in  the  Appendix  illustrates  this 
procedure  for  the  given  input  data  set.   Table  1 
displays  the  results  of  that  investigation  for  18  runs, 
including  the  first  already  mentioned.   Each  entry  is 
the  total  number  of  hour  colors  for  a  particular  run 
and  day  color. 

None  of  the  runs,  which  averaged  about  7  minutes  of 
IBM  360/67  time,  gives  the  acceptable  nine  or  less  hour 
colors  for  all  days.   Column  6  of  the  table  indicates  the 
number  of  acceptable  days  for  each  run.   Runs  3  and  18 
were  the  best  in  this  category  having  3  each.   Runs  1, 
15,  and  17  had  2  acceptable  days  each.   All  others  had 
one  or  less.   Column  7  of  the  table  indicates  the  number 
of  hour  colors  exceeding  9  each  day,  summed  over  the  five 
days  for  each  run. 
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4,   Conclusions. 

While  no  acceptable  coloring  solution  has  been  found 
directly  from  the  described  program,  several  program  runs 
are  sufficiently  close  to  indicate  that  modifications  in 
the  coloring  technique  may  in  fact  give  a  satisfactory 
solution  where  an  acceptable  trial  and  error  solution  is 
known  to  exist.   The  known  solution,  hand  generated  from 
the  same  data  with  all  constraints  taken  into  account, 
used  nine  periods  for  each  of  the  five  days.   Coloring  runs 
1  and  18,  Table  1,  both  came  within  3  periods  or  hour 
colors  of  equalling  the  number  of  periods  required  within 
one  week.   Runs  3*  12,  and  17  were  also  noticeably  close 
to  a  solution. 

In  this  case  where  a  solution  is  known  to  exist,  the 
results  obtained  must  be  attributed  to  the  coloring 
technique  applied.  With  an  acceptable  solution  in  mind, 
we  see  that  any  coloring  algorithm  or  technique  which  yields 
nine  or  less  hour  colors  per  day  is  a  satisfactory  pro- 
cedure.  Changes  in  the  present  program  in  DAYCOL  or  in 
the  algorithm  itself  could  improve  results.   One  such 
change  in  the  algorithm  may  be  the  following: 

T.  can  be  determined  by  the  Welsh  and  Powell  algorithm. 
The  set  of  nodes  in  T.  is  then  removed  from  the  graph, 
together  with  the  associated  edges.   The  remaining  sub- 
graph is  treated  as  a  new  graph,  and  the  Welsh  and  Powell 
is  reapplied  to  find  the  set  Tp ,  etc. 
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Also  a  change,  which  should  be  considered,  is  a 
coloring  program  where  the  day  and  hour  colors,  dependent 
upon  one  another,  are  assigned  together.   This  involves 
combining  the  functions  of  DAYCOL  and  HRCOLR,  with  what- 
ever algorithm  used,  into  a  program  operating  on  a  total 
conflict  graph  rather  than  separate  hourly  and  daily 
graphs. 

Another  observation  concerning  our  schedule  problem 
can  be  made.   Given  arbitrary  data,  we  may  question  the 
existence  of  any  schedule  solution.   A  solution  to  the 
presently  unsolved  problem  of  determining  the  chromatic 
number  for  an  arbitrary  graph  will  indicate  if  a  schedule 
solution  is  in  fact  possible. 
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APPENDIX 

GENERAL  COLORING  PROCEDURE 

FLOWCHART  A 


READ  IN  SECTION  VS.  COURSES 
AND  INSTRUCTOR  VS.  COURSES 
DATA 


DETERMINE 

CONFLICT 

LISTING 


ASSIGN  DAY  COLORS 
ARBITRARILY  USING 
DAYCOL 


READ  IN 
COURSE  HOURS 

MIA 


ASSIGN  HOUR  COLORS 
BY  ALGORITHM 
USING  HRCOLR 

PRINT  COLOR  VECTORS 
AND  NUMBER  OF  COLORS 
REQUIRE^ 
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READING  CONFLICT  DATA 
FLOWCHART  B 


DIMENSION  KX , INDEX , LIST , KD , NC , INDEXD , 
LISTD , IS ,NCOL , INDEXH , LISTH , R 
REAL*g  C,S,CC, BLANK, D,YKQ 
DATA  KX/Q/,INDEX/0/tLIST/0/ 


READ  (5,2) 
BLANK 
2  FORMAT 
(A7) 


DO  10  I  1,10000 
READ  (5,5)  K,S(I), 
(C(I,JJ  J  1,9) 


WRITE  (6,12)  IC1 
WRITE  (6.12)  IC3 
12  FORMAT  ( 5X,I5) 


WRITE  (6,13)  (S(I),(C(I,J), 

J  1  9),I«1.IC3) 

13  FORMAT  (5X.A5.q.9(lX.A7)) 
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COMPUTING  COURSE  LIST 
FLOWCHART  C 


CC(l)s 
C(l,l) 
N»l 


DO   20  1*1, IC3 
DO   19  J=1.9 


,|KX(I.J)»K 

7eS^KPsKP+1 


*KP*0 
NO  PO   18  Ksl.N 


-*i 


KN*N+1 
-*CC(KN)*C(I,J 
,KX(ItJ)*KN 


BVal 

N 

1 


19 


20 


VLC2 

N 


IP*13 

PO   100  Is l  IC2 
100  R(I)*RAN(IP) 


IC«IC2-U 

DO   110 
I»1.IC 


XMAXaR(I) 
JJ*I+1 
JKcI 


DO  108 
JaJJ, 
IC2 


XMAXsR(J) 
*JK*J 


.08 


XKQsR(I) 

R(I)=XMAX 

R(JK)SXKQ 


YKQsCC(I) 
CC(I)*CC(JK 

CC(JK)SYKQ 


DO  106  1*1, IC2 
DO  104  J*1.IC3 
DO  102  Ksl.9 
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PCXUKJsIU 


m 


06) 


WRITE  (6,12)  IC2 

WRITE  (6.16)  ((KX(I,J),J*1,9), 

16  FORMAT  (5X,6I9)/I-1,IC3 

WRITE  (6,32)  (I,CC(I), 
1=1, IC2)  32  FORMAT 
(5X,I5,2X,A7) 
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COMPUTING  CONFLICT  LISTING 
FLOWCHART  D 


— a— i 


|M4 DO  30  1*1,  IC2 

l| TINDKX(I]=M 


DO  30  Jsl, 
DO  28  Ksl.9 


IMsMrl 
INsINDEX(I) 


^ixM 


INDEX  (IC  2*1  )*M 

^FaM-1 

IC2PlaIC2*l 


DO  24  I  1,IC2P1,6 
11*1*5 

WRITE  (6,23)  (K,INDEX(K),Kc.I,II) 
23  FORMAT  (12(U,I4)) 


WRITE 

(6,12) 

M 


■M24 


bO  25  I  1,MF,6 

11*1+5 

WRITE  (6,23)  (K,LIST(K),KaI,II) 


■n25 
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REORDERING  CONFLICT  LIST 

FLOWCHART  E 


CALL  RELIST 
(LIST, INDEX, 
IC2) 


)0  3  IS1,IC2P1,6 

[P5*I+5 

*RITE  (6,23)  (K, INDEX ( 


K)>KSI>IP5)T^" 


DO  4  Isl,MF,6 

ip 5*1+5 

WRITE  (6,23)  (K,LIST(K),K»I,IP5) 


+i  4 


■■■)  - 


' 
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COMPUTING  UPPER  BOUND 
FLOWCHART  F 


INDEX  (IC2+1)*M 

DO  40  I»1.IC2 

40  KD(I)?INDEX(I+1)-INDEA(I) 


DO   51  Isl,IC2 
WRITE  (6,12) 
KD(I) 
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MAXVs 

0 

IC* 

Kfr-1 


DO  60  Irl.IC 

MAXsKD(I) 

JJsI+1 


DO  59 
J*JJ,IC2 


MAX*KD(J) 
JKsJ 


€H 


KQcKD(I) 
KD  ( I  )s  MAX 
D(JK)sKQ 
Klal 


MIN« 
KD(KI)+1 


WRITE  (6,64) 

MAXV 

64  FORMAT  (5X.I5) 
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READING  IN  COURSE  HOURS  AND  COMPUTING  D*Y  COLORS 

FLOWCHART  G 


DO  75  1=1,10000  READ  (5,65) 
K,D,(IS(J),Jrl  3)  65  FOR- 
MAT  (Ilt9XtA7>2X,3(lX,Il)) 


DO  73  M=l,3 

73  NC(M,L)SIS(M) 


©1 


WRITE  (6,66) 

D.(C(J,L),Jsl|9) 

66  FORMAT  5X.A7.2X. 3(11.11)) 


IC4TI 
*|WRITE  (6,12) 
IC4 


CALL  dAYcc!  (Nc, 

IC2.INDEXD, 
LISTD.NTOTC) 


■#»• 


WRITE   (6,23)  T- 

(K,INDEXD(K),Kslt6)| 


DO   SO  Isi.NTOTC.6 
IP5SI+5 
WRITE   (6,23) 
(K,LISTD(K),KgI^IP5) 


SO 
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COMPUTING  HOUR  COLORS 
FLOhGRAPH  H 


CALL  HRCOLR  (LIST, INDEX 
INDEXD , LISTD , INDEXH , 
LISTH.NCOL) 


WRITE  (6,6) 

(NCOL(I), isi,5) 

g  FORMAT  (5(51.15)) 


IWRltfi  (6,23) 

(KtINDBXH(K),Ksl.6) 


DO  g5  I=l,NTOTC,6 

IP5»I+5 

WRITE  (6,23) 

(KtLISTH(K),K*I,IP5) 
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FUNCTION  RAN (I) 
FLOWCHART  I 


FUNCTION  RAN(I) 
£1*65522 


TSTT 

2147463647 
*1 


RANsZ* 
.4656613E 

z2 


RETURN 
END 
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SUBROUTINE  RELIST 
FLOWCHART  J 


SUBROUTINE 

RELIST 

(LIST.INDEX.IC2) 


DIMENSION 

LIST 

INDEX 


DO  5 


INsINDEX(I*l)-l 
IM*INDEX(I)«\L 


DO  4  J*IK,IN 

MINxLIST(J-l) 

■'JKaJ-1 


KQsLIST(J-l) 

LIST(J-1)3MIN 
LIST(JK)aKO 


RETURN 
END 
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SUBROUTINE  DAYCOL 
FLOW GRAPH  K 


SUBROUTINE  DAYCOL 
(NC^Ca.INDEXD^LISTD^TOTC) 


DIMENSION  MA,ND,NC, 
JNDEXD,LISTD 
DATA  ND/0/tMA/0/ 
LDSl ___ 


IM*LD 
IN*LD*NC(1,I 


DO   1^  J?IM,IN 
LD*M0D(J-1,5)*1 
)-+ND(LD)?ND(LD)*l 
lMA(LD.ND(LD))SI 


DO   14  Kal.IN 

ND(LD)*ND(LD)+1 

MA(LD,ND(LD))*I 


DO  10  1*1.5 

INDEXD ( 1*1  )s  INDEXD ( I H-ND ( I ) 

DUND(I) 


DO   10  J31,IN 

LISTD(M)aMA(I,J) 

MsM-KL 


<S> 


NT0TC*M-1 

RETURN 

END 
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SUBROUTINE  HRCOLR 
FLOWCHART  L 


SUBROUTINE  HRCOLR  (INDEX, LIST, INDEXD,LISTD,NCOL, IN* 
DEXH.LISTH)  DIMENSION  INDEX, LIST, INDEXD,LISTD,INDEXH, 
LISTH^N^NDEXC^ISTC^DjKl^NZZ^COL^ISTTjNNZ^NNZZ 


DO  90  1*1,600 
INDEXC(I)*0 
90  KK(I)*Q 


DO  91  I»l,6000 
LISTC(I)»0 
91  LISTHap 


INDEIH(l)rl 
DO  30  m.5 


INDEXH ( 1*1 ):  INDEXD ( 1+1 ) 

IM*INDEX(I) 

IN»  INDEX  (MM 


M»0 
DO  29 
JSIM.IN 


Tdo 


0   2&  Kal.M 


INDEXC(K)sL 

IMl«INDEX(NN(K)) 

INlslNDEX(NN(K)+lK 


DO  26 
Il»l.» 


kpi?M^i 

INDEXC(MPl) 
:L 


VfHlB   (il32) 
I,M  32  FORMAT 
(//7X,I3.5X.I3)I 
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DO  22  J 
1,MP1,6 
JP57J+5 


WRITE    (6,23)    (K, 
INDEAC(K),KrJ.JP5) 
23  FORMAT    (12(4X,I4) 


22 


OJ 


Irl 


DO  20  J 
1,LM1,6 
JP5SJ»5 


fcRITE  (6.23) 
(K,LI3TC(K), 
KgJ,JP5) 


CD(J)sINDEXC(J*l) 
-INDEXC(J) 
CK(J)»J 


id1 


MMlrM-1 
DO  16 

Jsl.MMl 


WUsKD(J) 
->-JKrJ 
JJyJ+1 


eo_jhjk££ 


MAX»KD(K) 


H 


KQaKD(J) 

KD(J)«MAX 

KDjJKjsKQ 


KQ3KK(J) 

kk(j)skk(jk; 

KK(JK)5KQ 


DO  70   J*l,bOO 

NNZ(J)si 

DO  70  Kr 1.6000 


NZZ(J,K)r 


NZMlsNZ-1 

DO   31  KS1.NZM3 
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LlSTT(KK(J))sNCO 

NZZ(NZ)rJ 

NNZZ(NCO,NNZ(NCO))rKK(J) 


NZrNZ+1 
NNZ(NC0)*NNZ(NC0)«-1 
lMls lNDE^C (KK(J)) 


INialJ4L)EXC(KK(.)  j 

*1)-1 

JPl?J+l 


Juu  6  J1*JP1,M 
*HZMT»NZ-1 


-f 


2rlKDfiJt6(KK(Jl)) 
23lNDEIC(KK(Jl) 

Ifcl 


IN2=NNZ(NC0)-1 
DO  SO  K2SI2.N2 
DO  30  J2S1.IN2 


^    tel  LI5TT^X(Ji))sHC6 
(dOWNNZ(NZ)rJl 

V^y        N75N 


NZ^NZ+1 


NNZZ(NCO,NZZ 

(NC0))=KK(J1) 

11117.  ( [ICO )'  HNZ  ( NOP  )+l 


h/6Vnco^"^io)-> 


NC0(I)=NC0-1 

LLsi 

DO   50  JsIM.IN 


LISTH(J) 
LISTT(LL) 

LlsiL+1 


DO  UO  J*1.600 
INDSXC(J)sO 

KK(J)sQ 


1*9 


DO  60  Jsl, 

6000 

60  LISTC(J)=Q 


30 

RETURNI 
END 
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